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Mesoscopic non-equilibrium thermodynamics is used to formulate a model describing non- 

C/3 , homogeneous and non-Debye dielectric relaxation. The model is presented in terms of a Fokker- 

-)-* ■ 
2 ' Planck equation for the probability distribution of non-interacting polar molecules in contact with 

I ■ a heat bath and in the presence of an external time-dependent electric field. Memory effects are 

g : introduced in the Fokker-Planck description through integral relations containing memory kernels, 

I ^i ' which in turn are used to establish a connection with fractional Fokker-Planck descriptions. The 

model is developed in terms of the evolution equations for the first two moments of the distribution 

fT^ function. These equations are solved by following a perturbative method from which the expressions 

T-H ! 

t*~- ' for the complex susceptibilities are obtained as a functions of the frequency and the wave number. 

in 

1^ , Different memory kernels are considered and used to compare with experiments of dielectric relax- 

o 

j"*^ ' ation in glassy systems. For the case of Cole-Cole relaxation, we infer the distribution of relaxation 

times and its relation with an effective distribution of dipolar moments that can be attributed to 
different segmental motions of the polymer chains in a melt. 
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I. INTRODUCTION 

Relaxation processes in complex systems such as glass-forming liquids, polymeric melts, liquid crystals and even 
lipid membranes is a very active area of researcb^"— . They have been studied using different techniques ranging from 
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NMR and Brillouin light scattering to dielectric spectroscopj*ii^. The variety of relaxation phenomena occurring in 
these systems are non-Debye and therefore cannot be characterized by a single relaxation time'^''*. These multiscale 
time relaxations suggest an accumulative effect of different elemental relaxation processes that can be associated with 
collective motions in the system^Ti^. 

The average effects of these cooperative mechanisms observed in experiments can be accurately described by means 
of well known empirical relations such as the Williams- Watts function in time space, or through its conjugated 
frequency domain (Havrilak-Nagami) function. These functions, and many others that can be found in the litera- 
ture^rJl, reflect the specific features of the micro- or mesoscopic relaxation dynamics of each system. The microscopic 
description of these features may involve the formulation of specific molecule-molecule interaction models that lead to 
particular results. In view of this, it seems important to have a general theoretical model able to cope with different 
relaxation processes that may describe, for example, the dielectric relaxation in diverse materials including solids, 
glass- forming liquids and polymer melta^^. These last systems may exhibit, for example, anomalous Kerr effect which 
can be described in terms of models based on fractional Fokker-Planck equations (FFPE)ii^. 

In this work we adopt a mesoscopic point of view for describing the relaxation of a system of non-interacting polar 
molecules inmersed in a time-dependent external electric field, because at this level of description it is still possible to 
formulate the mentioned general framework able to cope with different relaxation dynamics. In particular, we use the 
mesoscopic non-equilibrium thermodynamics (MNET) formalism in order to derive a non-Markovian Fokker-Planck 
like equation (NMFPE) for the evolution of the probability distribution describing the orientation and position of the 
molecules. The MNET method uses the entropy produced by the system when it is driven out from equilibrium and 
is particularly useful for describing the stochastic evolution of systems at the mesoscopic scalei^"— . The advantage 
of this formalism lies in the fact that the NMFPE derived has a kinetic origin, and therefore there is no explicit 
assumption on the specific time correlation properties of the underlying Gaussian stochastic process, contrary to 
what happens with equations derived from Langevin formalism a^^'^" . Within this framework, linear relations between 
generalized forces and fluxes occurring in phase space are postulated, which allow us to introduce memory effects 
accounting for the non-instantaneous response of the system to the external perturbations. The effects of the non- 
Markovian dynamics on a measurable property of the system like its complex dielectric susceptibility can then be 
analyzed in a systematic way by considering different memory kernels. In this work we will concentrate mainly on 
memory kernels decaying as power laws. We will show that these particular functions turn out to be important since 
they lead to the mentioned FFPE. 



In the section II we will derive the NMFPE for a system of non-interacting polar molecules by using MNET. This 
equation describes both rotational diffusion and the diffusion of polarization, and we will show that it can be recast 
into a FFPE when proper memory kernels are considered. In section III we will construct the equations for the first 
two moments of distribution. Subsequently, in Section IV we will present a perturbation procedure which allow for 
a solution of the obtained dynamic equations, when the system is subjected to an external harmonic electric field 
superposed to a large static electric field. Consequently, we will obtain a closed expression for the non-homogeneous 
dielectric susceptibility of the system which incorporates non-Markovian effects, and the effects of the three first 
moments of the distribution function. These effects will be analyzed systematically in section V by considering 
diverse limiting cases and memory functions. Memory kernels decaying as power laws in time will be shown to 
generate FFPE and response functions of the Cole-Cole type^i. A comparison with recent experimental results is 
performed to validate the models derived. In section VI we will propose a model based on a theory of distribution 
of relaxation times which can be used to give an insight of the physical processes producing power law kernels and, 
accordingly, FFPE. Finally, we will present some final remarks an discuss the limitations of our model in section VII. 

II. NON-MARKOVIAN DESCRIPTION OF A SYSTEM OF POLAR MOLECULES 

We shall consider a system of N non-interacting polar molecules with dielectric moment, p, of magnitude p, which 
are in contact with a thermal bath at temperature T. For convenience, in the following we will use the unit vector, 
n = p/p, which represents the orientation of the molecular dielectric moment, and describe the system in terms of a 
single-molecule probability density. Let / (r, n, t) dfdn denote the probability for finding one molecule at a position 
ranging from r to r + df, and oriented within the range n and n + dn, at time t. 

The normalization condition, 

// dfdn f (r, n, t) = M, for ah t, (1) 

where M is the numerical density of molecules, imposes the following restriction on the time evolution of the probability 
density function fr^^ 

^ = -y^{m-R^{m^)■ (2) 



Here, Vi and Hi are the conjugate velocities to Xi and n^, respectively; Vi is the i-ih spatial derivative; and Ri is the 
i-th component of the so called rotational operator 



JXi — Sijk'IT'j 



d 

duk 



In this expression Sijk represents the Levi-Civitta antisymmetric symbol and summation over repeated indices will be 
implicit in the next. 

Consider now, that an external time-dependent and uniform electric field, E{t), acts on the the system. Then, the 
potential energy of one molecule in the presence of this field, U {n; t), is an explicit function of the orientation vector 
and also depends on time, t, since the electric field, E(t), is time-dependent. In first approximation, the potential 
energy of the system is given by 



U (n; t) = -pn ■ E (t) 



(3) 



which holds for molecules in which the anisotropy of the electric polarizability is negligible^i. 

The explicit form of Eq. ^ in terms of / (r , n, t) can be deduced by analyzing the entropy produced by the system 
during its time evolution. According to the so called Gibbs Entropy Postulat o^^i^^ , entropy is defined in terms of the 
probability distribution / by 



S = -kf 



drdn f (r, n, t) In 



/ {r, n, t) 



Cflcq 



(4) 



where kg is the Boltzmann constant and S'^'"^ is the local-equilibrium entropy. According to Eq. ([3]), the corresponding 
probability density in local-equilibrium is given by canonical distribution function 



n{n;t) 



zi°q (T; t) 



exp 



U{n;t) 
knT 



(5) 



where z'^^i is the corresponding partition function. 

In order to find the entropy production, a = d [S — S^"^) /dt, we take the time derivative of Eq. ([5]) from where. 



by using the normalization condition ([T]) and Eq. 1^, we find 



cr = l; ff drdn Lv, (/ V,) + nR, (f n,)] + kg II drdn / |- In /'°<i 



T 



mV. (/ V,) + Mi?. (/ ^.)J + kB J J drdn / - In /'°^ (6) 



where the non-equihbrium chemical potential, i-t — IJ. (r. n, t), has been defined by 

' / {r, fi, t) ' 



fj, (r, n,t) = ksTln 



/i-q(n;f) 



(7) 



Now, in accordance with the assumption of local equilibrium, we will consider external fields that evolve in a 
time-scale larger than the time-scale of fiuctuations. Therefore, changes of the local equilibrium distribution are 
not significant over the short time-scale of the fluctuations, and it can be approximated by f^°'^ {n;t) ~ /''^'^ {n;0). 
Consequently, the second term on the right hand side of Eq. ^ vanishes. This approximation is also supported by the 
fact that later on we will consider electric fields consisting of the superposition of a static, i?*^°^ , and a time varying, 
E^^' , contributions such that the magnitude of the time-dependent superposed field £^^'(t) is small when compared 
with the static field E^°\ that is, we wih assume E{t) = £:(°) +£^^'>{t) with |5(i)|/|i^^°^l << 1- 

Under this assumption, an integration by parts of Eq. ^ yields the following expression for a, 

a = -^ II drdn (/^^V,/x + fn,R,i?j , (8) 

where we have assumed an infinite system with vanishing probability density at the boundaries in r-space, and periodic 
conditions in n-space. 

The entropy production given by Eq. (|S]) involves the sum of products of the generalized currents, Vi and $1^, and 
forces, Vi/i and Ri^. Thus, following a linear response scheme^^, we may assume that the currents are proportional 
to forces and a linear relation for Vi and fi^, in terms of Vi/i and Ri^^ can be written. For simplicity, we will not 
consider cross terms in these relations, i.e. we will assume that the generalized force Vi/i is the unique responsible for 
producing the fiux Vi, and correspondingly for Ri^L and J7j. This assumption is in accordance with the corresponding 
Onsager-Casimir reciprocal relations given the vectorial character of the vectors involved in the description, see Refi^. 
In the general case and according to the phenomenology observed, we expect that the response of the system is not 
instantaneous, and then memory effects should be considered. Following RefSi ^^'^^ , these effects can be taken into 



account by introducing memory functions in the following linear relationships 



m 



ds K,J{t-s)f{s)yJ^l{s). 



(9) 



fQi = - ds Lij {t- s)f (s) Rjfj, (s) . 



(10) 



The memory functions Kij (t) and Lij (t), are such that Kij (t) = Lij (t) = 0, for i < 0. They describe the delay in 
the response of the system to forces which vary in time. This delay is always present and is particularly important 
when forces vary fast in time^^. 

By replacing Eqs. ([9]) and (ITOl) into Eq. ([2]), we obtain 



df /■' - /■* 

— = Vi / ds Lpij {t - s) \/jf (s) + Ri ds 4>ij {t - s) 



R^fis) + {^R,U 



(11) 



where for notation convenience, we have introduced the modified memory functions ipij (t) = ksTKij (t), and (j)ij (t) 
kBTLij (t). In order to obtain Ec^. dTTI) . we have also used the relations 



/V,/i = kBTVJ, and fR.^i = ksTRJ + fR,U, 



(12) 



which follow directly from Eq. ^. 

Eq. (jlip is the desired non-Markovian FP equation for the probability distribution function /. It provides a complete 
probabilistic description of the system at the mesoscopic level and constitutes a generalization of previous equations 
obtained in the literature for systems with Debye relaxation dynamics^^"^^. The memory kernels ipij (t) and (j>ij (i), 
characterize the translational and rotational diffusion processes in the system, respectively. For anisotropic molecules 
such as rigid polymers in solution, the translational diffusion tensor depends in general on the orientational degrees of 
freedom n^. Here, for simplicity, we will restrict ourselves to consider isotropic translational and rotational diffusion 
processes, i.e. (pij (t) — ip (t) Sij, and (f>ij (t) — (p (t) Sij. Consequently, our results will be expected to describe better 
the dynamics of slightly anisotropic polar molecules. Within this approximation Eq. (|lip can be recast in the final 



forr 



dt 



ds (p{t-s)f (s) + R, / ds cj){t- s) 



RJ {s) 



lis) 
knT 



R^U 



(13) 



In addition, notice that we have restricted ourselves to consider memory kernels which only depend on time, and 
accordingly, we have neglected the possible non-local effects from interactions at different points in space. 

In the following sections, Eq. (TTSl) will be used to obtain a dynamical model for the time evolution of the average 
polarization Pi and its associated order parameter, the quadrupole moment Sij. 



A. Homogeneous Approximation and Fractional Fokker-Planck Equation 

Diverse interesting limiting cases can be obtained from Eq. (J13p . Firstly, we will consider the so called homogeneous 
approximation?™, in which a simplified description in terms of the reduced orientational probability distribution, 



g (n, t) = df f (f, n, t) , 



is given. Averaging Eq. (|13p over space coordinates, the corresponding non-Markovian FP equation for g is found to 
be 



'^^=R.f_JsHt-s) 



kbT 



(14) 



Clearly, moments of the reduced distribution g will be uniform time-dependent quantities. Therefore, in the 
homogeneous approximation spatial diffusion of these moments is ignored and the properties obtained from Eq. (J14p 
must be considered as volume averages^^. 

It is interesting to notice that Eq. (|14p leads to a FFPE when the memory function (p (t) is conveniently chosen to 
be a power law. Specifically, let 



'{t) = { 



r(a 



'^{Dtf~^ fori>0 



(15) 



for t < 



where F represents the gamma function, a is a dimensionless parameter within the range < a < 1, and 13 is a 
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constant with units of inverse time that can be identified with the rotational diffusion coefficient of the molecules, 
as we will show below. It follows directly from Eq. ([M)) that when the memory function is given by Eq. ([Tsl) . the 
orientational distribution obeys the dynamic equation 



D" 



dg_ ^ 

dt ^ r (a - 1) 



Ri / ds 



1 



(t-s) 



2-Q 



knT 



(16) 



On the right hand side of the last expression we can identify the generalized fractional derivative operator- 



7.30 



.oJ^l'^g it) 



ds ■ 



9 is) 



rW 7 (t-s) 



1-1/ ' 



ioi v = a — 1. Thus Eq. (TlGl) can be recast in the form 



^ = D"^^l~'^R 



Rig 



ksT 



:R,.U 



(17) 



which is the referred FFPE for the orientational distribution. 

It can be noticed that Eq. ([T7)) reduces to the usual Fokker-Planck equation for a = 1 (normal rotational diffusion) , 



l^-A 



Rig + 



9 



knT 



RiU 



where D can be identified as the rotational diffusion coefficient^. For < a < 1, Eq. (|17p describes the non-Markovian 
behavior associated with rotational subdiffusive processes. 

An interesting point concerning the applicability of the memory function (|15p is worth stressing. The presence of 
r (a — 1) on the right hand side of Eq. ([T5|) yields a negative memory function for < a < 1. This implies that the 
memory function given by Eq. (jlSp superposes all the effects prior time t with an intrinsic negative sign. This fact 
is rather masked by the fractional derivative operator appearing in Eq. (J17p . and as far as we know it has not been 
noticed in literature where dynamic equations of the type (J17p are postulated for the probability distribution g. The 
presence of the function F (a — 1) is, however, necessary in order to obtain the FFPE, Eq. ([T6|) . from a power law 
decay of the memory function. 

Finally, we will consider the case in which both rotational and translational diffusion processes are present and can 



be described by power law decaying memory functions. For the translational kernel we introduce the function 



ifiit) 



r^l^i'3-2 fort>0 
for t < 



(18) 



where < /3 < 1, is a constant parameter, and Vp is a generalized translational diffusion coefficient with units of 
length square times t~^—. 

By substituting Eqs. ([T5]) and p^ into Eq. p^ . it can be easily shown that the FPE describing the evolution of 
the probability density / takes the form of a fractional equation, 



dl 
dt 



-.1-/3^72 



Vp.^'-^y'J + D^^^l-'^R, 



R^f 



f 



knT 



R^U 



(19) 



which includes anomalous diffusion through the fractional derivative -ooD^ . In the limit a,/3 = 1, the Markovian 
approximation is recovered. 

III. DYNAMIC EQUATIONS FOR THE MULTIPOLAR MOMENTS 

At the macroscopic level, the description of the system must be given in terms of the moments of the distribution 
function /. For this purpose, Eq. (jlip can be used to construct a hierarchy of equations for the multipolar moments 
associated with the charge distribution. Let us consider the dipole, Pi, quadrupole, Sij, and octupole, Qijk fields, 
which are defined by the following averages 



Piif,t) =JV dnf{f,n,t)p. 



(20) 



Sij (f, t) = -M^ I dnf{f, n, t) {SpiPj - p^Sij) , 



(21) 



Qijk {r, t) = -7V^ / dnf {f, n, t) bpiPjPk - 3p^ {Ip} 



fijk 



(22) 



In the last definition, {Ip} represents the symmetric part of the tensor Ip, which has the explicit form 
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{lp}»jfe = o (P^^jk + PjS^k + pAj) . 



Multiplying Eq. ([TT|) by Mpi and averaging over the orientational degrees of freedom we obtain 



m 

dt 



= V^ dsip{t-s)P^{s)-2 / ds <p {t - s) P^ (s) 



2Nf_ 
3k bT 



dscj){t~ s) El (s) 



SAfknT 



ds(f>(t — s) Sij (s) Ej (s) . 



(23) 



The first term at the right hand side of the last equation accounts for the diffusion of the polarization whereas the 
other terms are consequences of the orientational dynamics. 

Analogously, multiplying Eq. (jlip by Af^ [SpiPj — P^Sij) /2 and averaging over n, we obtain the dynamic equation 
for the quadrupolar moment. 



dt 



= V^ / ds ip [t — s) Sij (s) — 6 / ds(j){t — s) Sij (s) 



9Mp^ 
lOksT 



ds^{t~3) Pis)E{s) ~ / d3^{t-s)Q,jk{s)Ekis) 

LL Hi J 5A//crJ ./ 



(24) 



where P (t) E {t) is the symmetric traceless part of the dyad PE evaluated at time t, i.e. 



Pit)Eit) 



= E, (t) P, [t) + Ej (i) P, (t) - -Ek (t) Pk {t) 5,j. 

ij o 



It can be noticed from Eq. (P^ that the contraction Qijk (t) Ek (t) yields a symmetric traceless tensor of rank two. 
Consequently, Eq. (j24l) is in agreement with the symmetry properties of the order parameter tensor Sij . For notation's 
simplicity, in Eqs. (|23l) and (1^^ we have not written the dependence of the fields Pi, Sij and Qijk on r. 

Notice that Eqs. (|23p and ([M|) constitute the first two equations of an infinite recurrence hierarchy which includes the 
time evolution equation for the octupole moment, Qijk, and higher order moments. The problem of approximating this 
infinite hierarchy by a finite number of equations is a common feature of diverse systems in non-equilibrium statistical 
mechanics^i^^. In the following section we will describe in detail our closure approximation for the hierarchy of 
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recurrence equations for Pi , Sij , Qijk ■ 

In Eqs. ()23|) and (|24p . the moments Pt and S'ij turn out to be coupled only with moments of next superior order 
because translational diffusion has been assumed to be isotropic. Otherwise the evolution equation for Pi will not 
only involve Sij but Qijk as well, while the dynamic equation for Sij will include contributions of Qijk and the fourth 
order moment. These contributions have been shown to be important in recent literature for describing the dynamics 
of anisotropic molecules like those composing liquid crystals materials^^. 

Eqs. (l23l) and (p4)) can be reduced to the homogeneous case in which both Pi and Sij can be approximated by their 
position-independent volume averages 

P^{t)=M f dng{n,t)p„ (25) 



and 



S^j it) = ^AA2 f dng in, t) {ip,p, - pH,j) . (26) 



Then, it can be proved that P/' and ^Jl obey the dynamical equations 



t t 

- 2 



^^' ^2 / ds<^(t-s)P^(s) + |J% / ds(l)it-s)E,is) 



dt J ^ ' ^ ^ V / 3^^y 

— oo 

t 



iNkeT 



ds<t>it-s)SlAs)E,{s), (27) 



and 



t t 

^^" -6 ! ds4>it- s) Sf, is) + ^^^ I ds^it^s) \\P^ is) E (s) 



dt J ^ ' J 'i]\ J • lofcsT 



t 



bUkeT 



ds^it-s)Q,jkis)Ekis), (28) 



respectively. Eqs. ((27|) and (|28p can also be obtained from Eq. ([T4| for the reduced orientational probability distribu- 
tion after multiplying and averaging over n, as it was done in the case of Eqs. (|23p and (j24p . 
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IV. PERTURBATION THEORY AND DIELECTRIC RELAXATION 



In the following, we will focus on calculating the response function to a small oscillating field i?*^^-' (t) = ^(^^e*'^*, 
which is superposed to a large static field E^^\ such that £^^^ /E'-'^'f <C 1. Accordingly, the total field acting on the 
system is E (t) = E^^'' + f^^^e'"*. In this section we shall show that a perturbation theory in terms of the small 
parameter, £^^^ / E'^^\ can be used in order to find the frequency and field-dependent dielectric susceptibility of the 
system. 

With this purpose, let P^ , 5^- , etc., denote the electric multipolar moments associated with the static field. It 
is clear that these quantities relax towards equilibrium by following equations which are completely analogous to 
Eqs. dill) and dill), with pI^\ sf^\ Q^^l and eI°\ replacing the fields P„ S^j, Q^jk and E^, respectively. When 
both E^^^f and E'^^^ (i) act on the system, the response to the oscillating field can be identified from the dynamics of 
the perturbed fields P.^ = Pi — P^ , 5^- = Sij — 5^- , etc,^^. By substracting the equations corresponding to the 
unperturbed fields from Eqs. d^S}) and dMl), the following expressions are found for the time evolution of the perturbed 
fields 



ap 



(1) 



dt 



V=/..,(,-,)pf'(,„-2/<i,*(,-,„P.">(„ 



(29) 



-^ /-*"-')-.'"(») 



^ ' dscj,{t-s) S,,{s)Ej{s)-Sf{s)Ef 



'iNkeT 



dS. 



(1) 



dt 



= V 



t t 

2 ' dsip{t^s)S^^\s)^Q [ dscp{t-s)s\]\s) 



(30) 



IO/cbT 



ds (t>{t-s) \\P (s) E (s) - P(0) (s) i?(°) 



6 



dAfkeT 



ds<jy{t- s) \Q,^k {s) Ek is) - qI% (s) 4") 



For small perturbative fields, E'^^'' <^ E^'^\ deviations with respect to P^ , 5*1 , Q\jI, etc., are expected to be also 
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small. In this case, we can approximate 



S., (t) E, (t) S^ (t) Ef^ ^ 5(f (i) Ef [t) + S\f (t) Ef\ 



p (t) E (t) - p(") (t) ^(°'ll c^ [[p(o) (^) ^(1) (^)ii ^ rrp(i) (^) ^(0) 



fe it) E, it) - (M it) 4°) ^ q[% it) ii;W it) + Q^ it) 4°) 



ijk 



^ijk 



Moreover, for a large static field we can assume that the non-perturbed multipoles remain close to their equilibrium 



time-independent values, i.e. Pf> it) ^ P^ (e^'^A = P^, sf^^ it) ~ Sl^ f^(")) = S\ 



oq ^(0) 
■'ij ' ^ijk 



it) 



■'ijk 



(i?(o)) = 



*3i?fe' 'where, in order to simplify the succeeding discussion, we have omitted to write the dependence of the equilibrium 
moments on E^^' . In the following it should be regarded that these quantities are of order zero in the perturbative 
field. 

Finally, in order to close the hierarchy of recurrence equations, we will assume that moments of order Qijk and 
higher are not perturbed at all by the time dependent field, and consequently Q\^l it) ~ 0, for all t. 

All these approximations yield the following closed system of equations for the time evolution of the perturbed 



.(1) 



dipole, P- , and quadrupole, S*)- , moments 



(1) 



dP, 



(1) 



dt 



t t 

V^ f dsipit- s) PJ;^^ is) -2 f ds(l)it- s) P^'^^ is) 



(31) 



.^ /...(.-.)..-(.) 



iNkeT 



dscl^it^s)\s'r:tEf\s) + S^'\s)Ef 
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dS. 



(1) 



dt 



V" / ds¥.(i-s)4'^(s)-6 / ds(j>{t-s)si}\s) 



(32) 



lOfcsT 



ds(j){t — s) 



pcqj^il 



p(i)(s)J5(") 



6 



ds<P{t-s)Q-l,E\^\s). 



It can be observed that the sohition of Eqs. (I5T|) and ([5^ . gives indeed the perturbed polarization as a linear 
function of the perturbative field. From this solution the dielectric susceptibility of the system to the small oscillating 
field can be identified. For this purpose, we consider the Fourier transform of Eqs. ([5T|) and ([5^. which can be written 
in the form 



Gi 



{k,uj) pI^^ {k,uj) =2</)H|xo<5^ 



1 



3JVkBT 



S'^^E^ 



^)(fc,.)+<;)(fc,.)i^ 



(0) 



(33) 



Go 



{k,co)s<i^^{k,cu) = |xo0M{[[P'=^^W(fc,^)]]^^.+ [[pW(fc,^)^W 



6 



dMksT 



(/-H. 



(34) 



Here we have introduced the notation xo = Mp^ /iksT , for the static dielectric susceptibility, and Gi ( fc,w) and 
G2 ( fc, 1^ ) for the propagators in Fourier space defined by 



Gi 



I fc, w j = joj + 20 (w) + k^ip (w) , 



(35) 



G- 



I A;, oj ) = icj + 6(/) (cj) + k (p (uj) . 



(36) 



In Eqs. ([55]) and ([M]) we have also introduced the notation ii?*^^-' (A:, wj = (27r) £J'^^) (w)(5(fc), for the space-time 
Fourier transform of the perturbation field. 



The frequency-dependent dielectric susceptibility can be identified by solving Eq. (|34p for Slj , inserting the result 
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into Eq. ([33]) , and solving for P^ ' . After some algebraic manipulations this procedure finally yields 



Pi 



'^{k,Lo)^X^J{k,u;)E^'^{k,u:), (37) 



where the susceptibility, Xij {k,ujj, can be written as the product Xij — ^'hk ^kj, with the following definition for 



the matrices Mij and Nij, 



2 

M,, (fc, c.) = Gi (fc, uj) 5,, + ^^G^' {k, Lu) 02 (^) (3i?(°)2j,, + E^'^Ef) , (38) 



^B 



^4'^-) - ^(-M2^°^--wI;t^"+w4^^2"('''-)^^^)«^^^^"' 



„2 



^^ .G, 



'5k%T^^ 



'{k,u:)(t>H 



2 



^cq^(0)_^^cq^(0)^^_Z^oq^(„) 



(39) 



In the following section we will explicitely calculate the complex dielectric susceptibility, x \f^T<^) ^ ^or a specific 
geometry of the external fields and explore the effects of the memory functions, (uj) and if (oj), on this quantity. 

V. RESULTS 

In oder to illustrate the effects of the different mechanisms considered in previous sections on the response of the 
system to the applied fields, in this section we consider a specific geometry in which the static electric field, i?^°), is 
directed along the ea-axis of a Cartesian coordinate system. Furthermore, we will calculate the equilibrium tensors 
P°'^, 5,J"^ and Q'^'^f., by using the stationary solution of the reduced orientational equation, Eq. (ITi|) . 

p^(o) cxp(pn-^(")/fcBT) 
^ ^'^' ^ iTTkeT sinh(p^(o)/fcsT) ^^^' 

For this probability density and the geometry described above, we obtain the following expressions for the equilib- 



rium averages of the moments entering in Eqs. ([55]) and (|M| : 
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(PP)=pM 



( , \ 







(41) 



(5-)=/AA^(l- 



3i(0 



' -1/2 * 
-1/2 
1 



(42) 



(q:;^4°0-^'-^'^^°^ 



' -1/2 * 



..|^)mo-| 



V 



-1/2 
1 



(43) 



In Eqs. (|^T|) - (|43p . the dimensionless parameter ^ = pE^^^ /ksT quantifies the electric energy of a molecule in the 
static field relative to the thermal energy, and L (f) — coth(^) — 1/^ represents the Langevin function. As expected, 
from Eqs. ([55]). ([5^ and (|4ip - P5|) . we obtain that the matrices M,,,, Nij and Xij refiect the uniaxial symmetry of the 
problem and are given by 



' Afi ^ 



[M,, 



Mx_ 
M||y 



.{N,: 



'Nj_ ^ 
A^_i 

y iV|iy 



, and (xy) 



Xi 



After calculating the expressions for the components M±, M||, N± and A^|| (see the Appendix), the corresponding 
formulas for the components of the dielectric susceptibility perpendicular and parallel to the static field are found to 
be 



X-L 



(fc,w) 



1X0 



G^{k,u)+i^eG2'ik,u 



i(C) 



G^'{k,u;) 



e 



LiO 



(44) 






Gi(fc,w) + 1^2(5-1 (fc,^) 
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(45) 



where we have introduced the normahzed propagators Gi_2 — Gi,2 [k,uj] /20(aj). The dependence of Eqs. (|44l) - p5)) 
on the dimensionless parameter ^ imphes that the response of the system is in general a nonhnear function of the 
magnitude of the static component E^^\ This dependence also influences the characteristic relaxation time as known 
from experiments^. 

Our formalism reduces to well known phenomenological relations used to describe dielectric relaxation in diverse 
materials including solids, glass-forming liquids and polymer melt&i^, when the appropriate limiting cases are consid- 
ered. 

In particular, Eqs. pij) and psj) can be reduced to the Debye and Cole-Cole^i dielectric response functions if 
the infinite hierachy of dynamic equations for the moments of the distribution, starting with Eqs. (j29p and (j30l) . is 
truncated under the assumption that the order parameter tensor, Sij, remains very close to its equilibrium value, i.e. 



Sl'^' = S'^!^ ( E^^'' ] , and 5^- =0. If in addition we neglect the diffusion of polarization by assuming (p{t) = 0, for all 



t(0) _ QCq 

t, then Eqs. (|44l) and (|45|) . reduce to the simpler form 



" '"iw + 20((u;) 



where the amplitudes of the parallel and perpendicular components, Xn ^^'^ Xl; read^^ 



Xf = 3xo^p, and xl "= 2^^° ( ^ 



e 



(47) 



When the response of the system can be considered instantaneous, we have (f>{t) — DS{t), i.e. (f>{io) = D, and 
therefore Eq. (j46| takes the form of two independent Debye expressions 

Xi.||(^)=xl.||T-^^ , (48) 

where the relaxation time, tq, is related with the rotational diffusion coefficient through tq = 1/2D. On the other 
hand, for the power law memory kernel given by Eq. (|T5t . we have (p i^^) = -D" (*w) ", and Eq. ((46]) reduces to two 
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expressions of the Cole-Cole type 



X±,|| (^) = Xl.117— r Ta: 

■II 1 + (zWTcc) 



(49) 



where the relaxation time appearing in this expression, Tcc, is given by Tcc — l/2^/"Z3. 

In the following subsections we will analyze in more detail the dependence on the frequency, the normalized dipole 
electrostatic energy, the rotational diffusion and the diffusion of polarization of the real and imaginary parts of the 
complex susceptibility, x'±_ 11 = R-c {x±. || } and x" ,, = — Ini {x±. || }■ 

A. The case of Markovian dynamics 

For simplicity, in the following the analysis will be carried out in different successive reduced levels of description 
in order to better quantify and elucidate the effects due to couplings, diffusion of polarization and the non-Markovian 
character of the system. With this aim we have subdivided this section into subsections. 

1. Effects due to the coupling with higher order moments 

We will consider first the effects of the coupling of the average polarization with the moments of second and third 
order on the dielectric susceptibility. From Eqs. pSI) and ([55]). and (P^ - (|T7|) . it follows that coupling effects are 
manifested by the terms proportional to G2 ( fc, '^ ) , which is the propagator associated with the time and space 
evolution of Sij. Ignoring for convenience the effects associated to diffusion of polarization by making ip [lu) — 0, and 
the non-Markovian effects by assuming an instantaneous rotational diffusion kernel, i.e. <j) [uj) = D, the normalized 
complex susceptibilities perpendicular and parallel to £?(°\ Eqs. (HI)) and ps)) . are given by 



Gi{k,Q) =l + iG}, (50) 



G: 



(k,Oj) =3 + iw, (51) 
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10 10 10 

log ( (B ) 

FIG. 1: Imaginary part of the parallel component of the dielectric susceptibility for diverse values of the external field normalized 
with respect to the average thermal energy, ^. All curves correspond to Markovian processes in the homogeneous approximation. 
Solid lines contain the effect of the coupling of high order moments of the distribution function as obtained from Eq. (I45|l with 
Gi and G2 given by Eqs. (|5Up and (|5ip . Dashed lines are obtained from the imaginary part of Eqs. I|46p and l|47p . in which the 
coupling with high order moments is not taken into account (Debye case). 



where the normalized frequency, w, has been defined by 



w = Lj/2D. 



(52) 



When these expressions for Gi and G2 are used in Eqs. (j44|) and (j45|) . the result will contain exclusively the effects of 
the dynamics of the second and third order moments on x_l and xii. This result has to be compared with Eq. (j48p . in 
which these effects have been neglected. In Figs. [T] and [5] we present this comparison for the real and imaginary parts 
of the normalized susceptibility X|| (fc,wj — x\\ (^t'^] /xoi and four values of the normalized electrostatic energy, ^, 
ranging from small, S, — 0.5, to large values, ^ = 5, as compared with the average thermal energy per molecule. It can 
be observed that when the external static field is small, the coupling with high order moments produce no significant 
effect on the the complex susceptibility. However, for moderate, ^ = 1.0, and large external static fields, ^ = 5, the 
polar moments of the distribution increase and the changes induced by the coupling between these moments in both 
x' (oj) and x" ('^)i are appreciable. 

Figs. [Hand [5] show that the coupling of higher order moments reduces the susceptibility of the system with respect 
to the Debye-like case, in which this coupling is neglected. It can be also noticed that the frequency at which x'\\ takes 
its maximum value is shifted towards the region of large frequencies. However, the frequency dependence of x'\\ (^) is 
not modified in this limit since this quantity still decreases as ui~^ . 
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FIG. 2; The same as in Fig. [T]for the real part of the parallel component of the dielectric susceptibility. Solid lines represent 
the effect of the coupling of high order moments whereas the dashed lines correspond to the Debye case. 

2,. Effects due to diffusion of polarization 

The diffusion of polarization may also influence dependence on the wave vector and frequency of the dielectric 
susceptibility as we will show in this section. In order to focus on these effects, in this case we will consider the 
Markovian dynamics with memory functions 4> (w) = D, for the rotational diffusion, and (p (w) — P, for the diffusion 
of polarization, that has associated the diffusion coefficient V. We will also ignore the effects due to the coupling of 
the average dipole moment with moments of second and third order by making 6*2^^ = in Eqs. ((44]) and (|45p . 

Within these approximations, only the normalized propagator Gi is required for describing the behavior of X-LJIj 
and it can be written in the form 



Gi 



(fc,^) 



1 + "f + iuj, 



(53) 



where the parameter 7 — k^'D/2D, has been introduced to quantify the relative contribution of the diffusion of 
polarization with respect to rotational diffusion of the molecules. The normalized frequency uj is given as before by 
Eq. dSH). 

When Eq. ([5^ together with G2^ = 0, are used in Eqs. (|33]) and (|45l) . we obtain expressions which exclusively 
quantify the effects of the diffusion of polarization on the components of the dielectric susceptibility. Since in this 
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approximation the diffusion of polarization only modifies the real part of Gi , it is clear that its general effect on the 
imaginary part of x± a-nd x\\ "wiH consist in decreasing the value and shifting the location of the maxima of these 
functions in an amount proportional to the parameter 7. However, it can be noticed that the general dependence 
of the dielectric susceptibility with respect to the frequency is not modified by a memory kernel of the form ip — T), 
because both x'± and x\\ increase as w, for w <C 1 + 7, and decrease as tu^^ , ioi u ^ 1 + '^ . Similar conclusions are 
valid for the real part of x± and xii. 

B. The case of non-Markovian dynamics 

In this section, we will first consider the effects of the non-instantaneous response of the system to the imposed 
fields and afterwards both, the effects coming from coupling with higher order moments and those due to the diffusion 
of polarization. Consequently, we set (f (uj) = 0, and Gj"^ = 0, in Eqs. pi)) and P5|) . This yields Eqs. (05]) and (|T7)) 
for the dielectric susceptibility components, in which the memory kernel, '/>(a;), has to be given in order to specify 
the behavior of the dielectric relaxation. 

For the sake of clarity, in the present paper, we will only analyze the effects of two types of memory kernels, namely, 
the power law memory function given by Eq. (J15l) . and a modified kernel that combines both exponential and power 
law decay. We will write the latter in the form 

^.^e-2^* (2^i)""' fort>0 
'(t)^{^^-'' , (54) 

for t < 

with < a < 1. 

Let us discuss first the effects of the memory function given by Eq. (1151) . Previously, we have shown that this 
function yields the FFPE, Ec^. p7)) . and the Cole-Cole law for the dielectric susceptibility, Eq. P^ . Accordingly, the 
effects of this kind of non-Markovian dynamics can be analyzed by simply comparing our previous results, namely, 
Eq. ((33) with respect to the Debye expression, Eq. (|T7)) . in which such effects do not appear. 

In Fig.|3]we perform this comparison for the real and imaginary parts of x\\ i^)- For convenience, we have illustrated 
our results in terms of the normalized susceptibility x\\ = X\\/Xo, and the normalized frequency iD, Ec^. (1521) . 

Curves in Fig. [3] exhibit the well known features of the Cole-Cole susceptibility. The dielectric peak in the Cole- 
Cole case is symmetric and behaves in the form x" ^ "^"i x" '^ '^~"i in the limits of small and large frequencies, 
respectively. For < a < 1, the width of the dielectric loss increases the smaller a is, while the departure with respect 
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FIG. 3: Normalized parallel component of the dielectric susceptibility in the non-Markovian homogeneous approximation. The 
memory kernel describing the rotational diffusive processes is considered to decay according to the power law Eq. (|15|) . The 
dimensionless exponent a, takes diverse values including q = 1, which corresponds to normal diffusion, and a = 0.9, 0.8, 0.7, 0.6, 
which describe rotational subdiffusive processes. Xii' S'^d Xii ^re calculated from Eqs. (|49[) and (|48p . with X|| = Xjl/Xo- The 
value of the normalized field is taken as ^ = 1. 



to the Debye case is eliminated for a = 1. The frequency dependence of the real part of the susceptibility is also 
modified by the power law memory function, Eq. (jlSp . and instead of the characteristic Debye's decay, %' ^ w^^, a 
decay of the form x' ^ w^", is observed. 

The modified memory function given by Eq. (|54l) yields an asymmetric shape for the dielectric peak. This result 
follows from the Fourier transform of Eq. ()54|) : (w) = D (1 + iu)/2D) , and after replacing this expression into 
Eq. (|46p . The resulting formula for the normal and parallel components of the susceptibility can be written in the 
form 



X-L.ll (^) = 



^III 



1 + IUJTq (1 + iUJToY 



(55) 



where we used tq = 1/2_D and the amplitudes x°l ii ^'"S given by Eq. (|47| . 

In order to analyze the low and high frequency behavior of Eq. ([55]) , we will use the formula 



(1 + iwTo)"^' = (l+wV2)^e*("-i)*"""'(""°). 
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log (CO) 

FIG. 4: Normalized parallel component of the dielectric susceptibility in the non-Markovian homogeneous approximation. The 
memory kernel describing the rotational difTusive processes is considered to decay according to the modified power law Eq. (|54[). 
For a — 1, the usual decay in the dielectric loss Xii' ~ w , is observed. For a < 1, x" becomes asymmetric, increasing as to 
and decreasing as (jj~" . x\\ a-nd x'\\ ''■r^ calculated from Eqs. (|55|l and (|48|) . with x\\ = X||/Xo- The value of the normalized field 
is fixed at ^ = 1. 



In the limit lutq <C 1, the imaginary part of X±,|| increases linearly with w, i.e. X±.\\ ~ '-^'i while for ujtq ^ 1, it 
decreases as X±,\\ ^ w^". This asymmetry in the dielectric loss is not exhibited by the Debye nor by the Cole-Cole 
functions, but resembles the dependence that is observed in systems described by the phenomenological Cole-Davidson 
expression!'!. In a similar way, it can be shown that the real part of Xj_.|| decreases as X±.|| '~ w~" for lutq ^ 1, in 
contrast with the usual Debye dependence. Fig. |4] exhibits these effects on the real and imaginary parts of x\\ ('^); in 
terms of the normalized susceptibility X|| — X||/xoj and the normalized frequency w, Eq. (j52p . 

The general case in which both, rotational diffusion and the diffusion of polarization processes are present and 
have a non-Markovian character can be analyzed as a function of the parameter ^. In general terms, the qualitative 
correction arising from the coupling with higher orders moments only modulates the amplitudes of the obtained 
functions and their characteristic relaxation times (their location along the w-axis), as shown in Figs. [T] and [2] Thus, 
if for simplicity we consider small values of ^, such that the coupling of high order moments is neglected, then the 
expression for the components of x {k,io] takes the form 



^^'" ^'^^ - ^^■llzc. + 20(c.) + fcV(c.)' 



(56) 
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with x']_ II given by Eq. (j47|) . When the memory kernels (f) (t) and ip (t) are assumed to decay as power laws in time, 
Eqs. (US]) and ^, then Eq. ^ reduces to 

^ ^ 2 \ D ) ^ 2 VDJ 

which, as Eq. (|55l) , is also asymmetric. Unlike in the previous case, now the low and high frequency behaviors of 
the real and imaginary parts of the susceptibilities will be determined by the exponents a and /3. This result implies 
that the diffusion of polarization introduces asymmetries in the response functions of the system that depend on the 
viscoelastic properties of the medium. The representation of Eq. (|57p contains a combination of the effects shown in 
Figs. [3] and H 

C. Cooperative effects 

A simple and effective generalization of the present model can be formulated along the following lines. Consider the 
case of a system in which different kinds of polar elements, characterized by their dipolar moments, p^'^' , are present. 
In this case, the probability distribution function will depend on several orientation vectors, n^'^' , associated to each 
type of polar element. For example, in the case <z = 2, the distribution function f {r^^\ r^^\ n^^\ n^^\ t) , will satisfy 
a bivariate conservation equation entirely similar to Eq. ([2]), but having the corresponding divergence term in the 
^^1) ^ ^2) ^ j^(i) ^ ri(^))-space. 

As a consequence, a bivariate Fokker-Planck equation similar to Eq. (jl3|) can be obtained from the mesoscopic 
thermodynamics analysis. This equation can in turn be used to formulate the corresponding evolution equations for 
the moments of the distribution, in which P^' ij^t) = Mq J dn'-^^dn'^^'>dr^^'>dr^^'>pi'' fS{'r~r^'^'>), q — 1, 2, are the dipole 
moments, see Eq. (l20l) . Here, J\fq is the number of polar elements of class q, and the operator (5(r — f^'''), evaluates 
the value of the field P^"^' in the continuous position r, see for example^. The definitions of the quadrupolar and 
octupolar moments are therefore given by entirely similar relations as those of Eqs. (|2ip and (j22p . 

If we assume that these polar elements are independent, a simple and powerful model can be formulated to account 
for experiments. Under the independence assumption, the evolution equations associated with P''^' (r , t) and P^^' (r , t) , 
are also independent and therefore, the total relaxation process results from the sum of both contributions. Taking 
into account only the time dependent field, the relation between the polarization and the applied electric field can be 
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written as 



P^^i2 ^i"^(^''^) = Xv{Kuj)E,{k,uj), (58) 

ni—l 

where Xij{k,oj) — '^m=iXiT {k,oj), since in the present approximation the response function is independent of the 
apphed field. 

As an example, consider that q = 2 and that the diffusion of polarization effect is subdiffusive anomalous. Thus, 
by applying a similar procedure to the one yielding Eq. (j57|) . we obtain 

where the subscripts for the parallel and perpendicular components have been suppressed in order to simplify the 
notation. In this case, the susceptibility is also asymmetric but it may present two broad peaks when the characteristic 
relaxation times associated to P*-^' and P^^^ are sufficiently separated. 

The imaginary part of Eq. ([55]) can be used to account for the dielectric relaxation spectra of glassy materials. 
This is shown in Fig. [5l where a comparison between experiments (symbols) and theory (lines) is presented in a 
normalized representation of the dielectric loss, x", see Refi^. The experiments were performed with 50/50 (wt/wt) 
glycerol/water samples at 160 K (solid), in the absence and the presence of myoglobin at different levels of hydration. 
The symbol h stands for the level of hydration, where h = oo corresponds to the case when no myoglobin was added. 
The dielectric relaxation spectroscopy experiments were performed by introducing the samples in a capacitor and 
applying an oscillatory electric field''®. The presence of the hydrated myoglobin broadens the spectrum and shifts 
the maximum towards low frequencies. In Fig. [5] we have included two cases: (a) Without considering diffusion of 
polarization, i.e. V^-^ = Vj^^ = 0, in Eq. (|59|) . (b) Considering the diffusion of polarization, i.e. P^j ^ V^^ ¥" 0- It can 
be observed that the diffusion of polarization improves the agreement between theoretical and experimental results. 

It is interesting to notice that the diffusion of polarization appears to have most dramatic effects in the absence of 
myoglobin, while its effects are small in the presence of myoglobin. These effects can be understood by considering 
that the reduction of the amplitude of the dielectric loss is due to the presence of both the myoglobin and the hydration 
shell, whereas the addition of new modes at low and high frequencies can be attributed to the presence of the hydration 
shelP®. Then, it is plausible to assume that the hydration shell induces long-ranged orientational correlations between 
solvent particles, which reduce the relative importance of the diffusion of polarization term with respect to the non- 
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FIG. 5: Nonnalized dielectric loss of glycerol/water samples at 160K in presence of myoglobin with different levels of hydration. 
The symbols represent experimental results taken from Ref.^^ whereas the lines are the fits from Eq. ()59[) . Our results 
suggest that the structural a-relaxation process affects the observed /3-relaxation spectrum by introducing an asymmetry in 
the amplitudes of the spectrum at low frequencies. 

diffusive terms in all the evolution equations for the multipolar moments. The orientational correlations explain why 
in Fig. [5ji the agreement for the larger value of the hydration number (h — 2.5) is better at low frequencies than 
in the case of low hydration numbers {h ~ 0.5). For low h, the hydration shell is less important and therefore the 
amplitude of the dielectric loss increases as well as the importance of the diffusion term in the evolution equations for 
the multipoles. 

In this specific application, the superposition of two independent dielectric responses given by Eq. (|59p , accounts 
very well for the complicated shape of the dielectric loss found experimentally. According to our model, this function 
should consist of a peak centered at small frequencies, associated with the so called a-relaxation processes, superposed 
to a peak located at larger frequencies, which is associated with slower or /3-relaxation processes. The latter peak 
is the one observed in Figs. [5] (a) and (b). It can be noticed that the diffusion of polarization allow for fitting 
the experimental data in the range of small frequencies. Notice that, as h increases, the exponent /? increases as 
well, implying that 13(h) and indicating that hydration favors a transition from anomalous to normal diffusion of 
polarization, see the tables in Appendix B. 

VI. DISCUSSION 



A physical insight of the mechanism underlying the power-law behavior of the memory kernel in Eq. (jl5|) can be 
elucidated by decomposing the corresponding spectrum in terms of elementary Debye processes, that is, by appealing 
to the theory of distribution of relaxation times^. For simplicity, our analysis will be focused on the description 
leading to Cole-Cole parallel susceptibility (PU|) with the use of the power law memory function, Eq. ([T| 
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We will proceed by assuming the existence of a distribution of elementary processes with proper relaxation times 
and amplitudes, whose superposition generates a wider spectrum such that corresponding to the Cole-Cole function, 
Eq. (j49[) . In this form, the Cole-Cole function can be recovered from a superposition of elementary Debye processes 
through the relation^ 

XII H= r d(ln\f\),p(ln\f\)—^= r dy^Piy)-—^ -, (60) 

where for simplicity in notation we have introduced the variable f = t/tcc and used the definition y = In [f| in the 
second equality. The function ■0 (y) is a continuous distribution of relaxation times having the explicit form 

X|| sin (an) 

27r cosh (ay) + cos (an) 



In fact, the analysis can be simplified by considering a discrete expansion of Eqs. (|60| and (|6T|) . that is, by 
approximating the Cole-Cole response law by a superposition of independent Debye-like processes Xk (w), distributed 
over a discrete space in the form 

OO OO J 

XII H= Y. »H= E T-^' (62) 

fc— — OO fc— — OO 

where the amplitudes, Ak, and relaxation times, Tk, of the fc-th process are given, respectively, by 

A _ '''^11 sin(a7r) j _ fc ^co^ 

27r cosh (afc) + cos (aTr) ' 

Notice that the index k in the summation of Eq. (1631) , quantifies the time-scale separation between the fc-th elementary 
process and the process with the relaxation time Tcc, which defines the location of the maximum of the dielectric 
response function. Those processes labeled by an index fc 3> 1, evolve toward equilibrium in a characteristic time 
which is much larger than Tcc, while those with fc <C — 1, decay much faster than the processes with relaxation time 
of the order of Tec- 
Fig. [S] illustrates the approximation of both the real and imaginary parts of the normalized Cole-Cole function. 



X|| i^) _ 1 + (tccUj)"' cos (a7r/2) 



Xy [1 -I- (tccw)"] -I- {tccUj) " sin^ (a7r/2) 



(64) 
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FIG. 6: Approximation of a Cole-Cole function by a superposition of discrete Debye components with amplitudes and relaxation 
times given by Eq. ((63}. The real and imaginary part of the Cole-Cole function, Eq. (|49|) . are represented by the continuous 
curve. They correspond to a = 0.6. The superposition of Debey components, Eqs. (|62p and H63p . is represented by the square 
symbols and it can be observed that it fits the Cole-Cole function over the whole range of frequencies. Some of the Debye 
components, those with indexes k = —2, —1, 0, 1, 2, are explicitly shown. 



X|| (^) _ {TccUj)°'ain{aTT/2) 

X|j " [1 + {T^^ujrf + ir^^Ljf" sin2 {aTT/2) ' 



(65) 



by the superposition of Debye components over the discrete distribution of relaxation times given by Eq. (|63p with 
a = 0.6, and vifhere the infinite sum in Eq. (j63p has been approximated by a sum over a finite interval ranging 
from k = —100 to k = 100. Eq. ([M]) constitutes a very good approximation to the Cole-Cole function, reproducing 
very well the behavior of the latter over the whole frequency range considered. It is interesting to notice that a 
partial superposition of independent Debye contributions in which elementary processes within a given range are 
suppressed, yields a two step relaxation dynamics manifested through the dielectric loss with well marked peaks 
as those observed experimentally"^^. This elimination of elementary processes can be represented by non- uniform 
distributions of characteristic times^. 

These results can be used to give an insight of the physical mechanisms underlying the relaxation process if we 
consider Eqs. (|Tfl) and (05]), and the definition of the static susceptibility xo — Mp^ /ikBT. The amplitude of each 
elementary Debye contribution in Eq. (|62]) , can be written in terms of an effective numerical density of polar elements 
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A/fe participating in the process with relaxation time r^ , that is 

ksT 5 

Comparing Eqs. (I63p and (|66p we obtain 

^ ^ sin {an) 

6ttL (f ) cosh (ak) + cos (an) ' 

where Uk = Afk/M, is the effective fraction of polar elements which contribute to process k. In the limit of slow and 
fast processes, i.e. A; ^ 1 and /c <C — 1, respectively, and in terms of the relaxation time r^, Eq. ([57]) . n^ can be recast 

into 

£ f;^) for Tfc < Tec 



This result implies that the power law memory kernel (ITS]) and consequently the Cole-Cole function P^ . can be 
suitable for describing the dynamics of systems in which a distribution of relaxation times exist, such that the fraction 
of polar elements of the system participating in the processes with relaxation time r/j, scales as r^ in the range of fast 
processes, while in the range of slow processes it scales as TjT". In addition, it should be noticed that the fraction of 
polar elements participating in a given relaxation process depends on the dimensionless parameter ^ which in turn 
depends on the static applied field E^^\ 

VII. CONCLUSIONS 

Starting from purely thermodynamic grounds, we have formulated a very general framework, which is useful for 
describing non-Debye dielectric relaxation in diverse materials including solids, glass-forming liquids and polymer 
melts. The framework is based on the derivation of a non-Markovian Fokker-Planck like equation that incorporates 
memory effects for both rotational diffusion and diffusion of polarization and therefore provides a very wide range 
of possible behaviors as analyzed through section V. We have shown that fractional Fokker-Planck descriptions can 
be derived as particular cases of our model when a specific form of the memory kernel is chosen. The formulation of 
the model on thermodynamic grounds simplifies its generalization to cases in which external temperature or pressure 
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gradients, among others, are applied on the system. 

The evolution equations for the polarization vector and higher order multipoles can be derived systematically. 
We analyzed the evolution equations for the first two multipoles by closing the corresponding hierarchy assuming 
that the higher order multipoles may be considered near to equilibrium. This assumption yields previously known 
results for the dependence of the susceptibilities and the relaxation times when the system evolves in the presence of 
a time dependent externally applied electric field consisting of a static and a harmonic dependent component. We 
analyzed in detail the effect that each mechanism considered has in the behavior of the susceptibilities and shown how 
empirical formulas can be obtained from the general formalism. In addition, we outlined how the present approach 
can be generalized to consider the non-Markovian dynamics of several polar elements and the effect of this on the 
total susceptibility. In this case, our results were used to fit with recent experiments performed with glassy materials 
damped with myoglobin proteins at different levels of hydration. The agreement between theory and experiment was 
very good. 

The analysis of the susceptibilities by means of the theory of distribution of times lead us to propose a physical 
interpretation of each relaxation process in terms of the fraction of polar elements that participate in the relaxation 
process for a given frequency of the external field. We found that this fraction of polar elements, related to a 
particular characteristic relaxation time, depends on the static externally applied field E^^' . This result indicates 
that the application of a static field may be used to control the fraction of polar elements that participate in a given 
relaxation process (possibly including multistep relaxation) and thus can be useful in both the analysis and the design 
of materials with specific absorption windows. 
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Appendix A: Auxiliar formulae 



In section V, we have discussed how the uniaxial symmetry of the problem reduces the matrices Mij , Nij and Xij to 
diagonal forms that contain two contributions, one parallel and one perpendicular to the static applied field. Making 
use of Eqs. ([55)l . (|5^ and (|lT |l - ([i5)) . it can be shown that the explicit form of the elements M±, My, N± and TVy are 
given by the expressions: 



M^(k,oj]^Gi(k,Lj]+-C 



G2(k,io 



(Al) 



Ml 



ll(fc,c.)=Gi(fc,c.)+^C 



5 G: 






(A2) 



-(-)--(H^M-^)-ig[-(f^i) 




LiO 



(A3) 



iV|| i^k,ujj = 6x0(1} [k, I 



Mi) 
^ ■ G: 



2(j)(k,Lu\ 



-1 + ^^(0 



(A4) 



From these equations, the explicit form of the susceptibilities Eqs. p4)l and ((45)) was calculated 
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Appendix B: Fitting parameters 

The parameters used for fitting the experimental data with Eq. (j59p and represented in Figs. [S] (a) and (b), are 
jivcn in the following table. 



Model 


X° 


X? 


Ql 


02 


/3i 


02 


Tl 


1-2 


^/3l 


■•■02 


Fig. [11(a) 


1.368 


0.859 


0.715 


0.374 






0.063 


94.470 






Fig. [11(b) 


4.408 


0.978 


0.831 


0.013 


0.703 


0.541 


0.177 


0.673 


0.707 


0.421 



We have introduced the notation Tq = l/2Dq, rp^ — l/k'^Vp^, for g = 1, 2. 



